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Abstract 

The dynamics for a thin, closed loop inextensible Euler's elastica moving in three 
dimensions are considered. The equations of motion for the elastica include a wave 
equation involving fourth order spatial derivatives and a second order elliptic equa- 
tion for its tension. A Hasimoto transformation is used to rewrite the equations in 
convenient coordinates for the space and time derivatives of the tangent vector. A 
feature of this elastica is that it exhibits time-dependent monodromy. A frame frame 
parallel-transported along the elastica is in general only quasi-periodic, resulting in 
time-dependent boundary conditions for the coordinates. This complication is ad- 
dressed by a gauge transformation, after which a contraction mapping argument can 
be applied. Local existence and uniqueness of elastica solutions are established for 
initial data in suitable Sobolev spaces. 
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1 Introduction 



1.1 Statement of the result 

Let x(s, t) be a smooth closed curve in R 3 , parametrized by its arc length < s < 2n 
and by time t. Physically, the curve can be thought of as a loop of very thin inextensible 
wire which can move in space. We assume that the curve is flexible, with potential energy 
at time t determined by its curvature k(s, t) = \d 2 x(s, t)\, 

i r 2w rh 
V = 5 / K 2 (s,t)^. (i.l) 

Its kinetic energy at time t is given by 

r = \fj\ d ^)\ 2 ^- 0-2) 

We will refer to a dynamical curve x(s, t) with this potential and kinetic energy as a 
dynamical elastica. This choice of kinetic and potential energy gives rise to the variational 
problem of finding extreme solution curves for the Lagrangian C = T — V subject to 
the constraint that s be in fact arc length. The constraint can be implemented by adding 
a suitable Lagrange multiplier term to the Lagrangian, enforcing the condition that the 
tangent vector u(s, t) = <9 s x(s, t) be a vector of fixed modulus one. The variational 
problem we consider is then 

5^|£(x(-,t),^x(-,t))-i^A( S ,t)(|a s x( s ,t)| 2 -l)^} = 0. (1.3) 

The variational derivative with respect to A provides the arclength constraint. The result- 
ing Euler-Lagrange equations are 

d 2 x = -dt* + d s (\d s x) , |<9 s x| 2 = 1, (1.4) 

where, we emphasize, the Lagrange multiplier A(s, t) is a function of both s and t. From 
this equation, it is apparent that A(s, t) has the physical significance of being the tension 



of the curve at s, at time t. Note that A can be negative, signifying compression. Eq. ( |1 .4[ ) 
is nonlinear in x due to the fact that at any time t, the tension is itself the solution to 
an inhomogeneous elliptic equation involving the tangent vector u(s, t) = d s x(s, t) and 
its derivatives (see Eq. ( |1.21| ) below). We note that omitting or relaxing the arclength 
constraint would require computing the curvature as k = d s x x d 2 x/\d s x\ 3 , leading to a 
rather forbidding quasilinear Euler-Lagrange equation. 

Differentiating Eq. Q1.4D with respect to s, one obtains for the tangent vector the equa- 
tion 

d 2 u = d 2 s (~d 2 s + A) u , u G § 2 . (1.5) 
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We study here the Cauchy problem for the wave equation Eq. (|1.5|), assuming periodic 
boundary conditions 

u(s + 27r,t) = u(s,t) , (1.6) 

and initial conditions 

u(s,0) = uoO) , dt\i(s, 0) = ui(s) , (1.7) 

where u and and Ui are given 2tt -periodic functions with 

|uo(s)| = l, uo(s) • Ms) = (1.8) 

for all s, that is, u (s) lies in § 2 , and Ux(s) is a tangent vector to the sphere at the point 
u (s). We are interested in weak solutions of Eq. (|T3|), defined by the property that 

-jf </>(*) -<9,u(M)^ = y 9 2 0( S )-(-9 2 u( S ,t) + A( S ,t)u( S ,t))- (1.9) 

for any smooth 27r-periodic function (p. We require the map t t— > (u, 9 t u) to be strongly 
continuous in a suitable Sobolev space, and interpret the time derivative on the left hand 
side in the sense of distributions (see Section for a definition of the Sobolev spaces). 
The following is our main result. 

Theorem 1.1 Let (uo, ui) be a pair of 2ir -periodic functions in a Sobolev space H r+2 x 
H r with values in M 3 which satisfy the constraint in Eq. ( |i.#p . For r > 1/2, there exists 
a time T > 0, which depends on ||u || Hr+ 2 and on HuiH^, such that the initial-value 
problem given by Eqs. (|0|)-(|77^) has a strongly continuous solution u on [0, T] with 
u(-,£) G H r+2 , dtu(-,t) G H r , and X(-,t) G H r+1 . The solution is unique and depends 
continuously on the initial data. The conclusions hold for allr > if the initial values Uq 
and ui lie in a common plane through the origin. In this case, the solution is planar. 

Theorem |1 . 1| implies that Eq. ( |T~4t ) with initial values 

x(s,0) = Xo(s), $x(s, 0) = Xl (s) (1.10) 

satisfying the compatibility condition 

<9 s x • <9 s xi = , (1.11) 

is well-posed in H 7 ^ 2 x if 3 / 2 (in H 3 x if 1 , if the motion is planar). Periodicity of (x , Xi) 
implies that x(s + 27T, t) — x(s, t) and its time derivative are zero at t = 0, whereas 

/■27T 

a 2 (x(s + 27r,t) -x(s,t)) = / 9 2 u(s',t)ds' = (1.12) 
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by Eq. (JT3|), so that the curve x(-, t) remains a closed loop. 

Theorem [TTT] applies also to infinite (expanding or contracting) helical curves x(s, t) 
for which u(-, t) is still periodic. Eq. ( |1.12[ ) implies that unless ui = d t u(-, 0) averages to 
zero, x(s + 2tt, t) — x(s, t) grows linearly in time. Since it cannot grow beyond 2n, the 
curve must disintegrate in a finite time under these circumstances. 

Since the total energy T + V is equivalent to the natural norm of (x, <9 t x) in H 2 x L 2 , 
it seems reasonable to consider the Eq. fll.4| ) in H 2 x L 2 . By conservation of energy, a 
small-time existence result for solutions of Eq. ( |1 .4[ ) in H 1 x L 2 would imply that solutions 
exist globally in time. This would amount to proving Theorem |1.1| with r = — 1 for initial 
values (u , ui) where ui averages to zero. It is an open question under what conditions 
on the initial values the solutions to Eq. ( |L4] ) exist globally in time. Blowup, if there is 
any, must involve a transfer of energy to the high-frequency (large n) modes of the Fourier 
transform of the solution. 

Our theorem should be compared with results of Caflisch and Maddocks [(TJ], who 
have given a proof of global existence for a planar dynamical elastica. Their equations 
of motion include an additional rotational inertia term proportional to | J \d t u\ 2 ds in the 
kinetic energy T (see their Eq. (2.16)). They assume that the initial values u and ui 
are piecewise C 2 and piecewise C 1 , respectively, a stronger assumption than we require 
in the planar case (although it appears that their smoothness assumptions can be relaxed 
somewhat). With this additional kinetic energy term, the equation of motion (for the 
tangent vector angle) can be written as a single second order semilinear wave equation 
with a non-local nonlinear term. Conservation of energy ultimately provides the a priori 
bounds needed to prove their global existence result. 

1.2 Related geometrical and physical problems 

The variational problem stated in Eq. Q1.3D is closely related with wave maps. For a wave 
map problem, a typical choice for the kinetic energy would be 




(1.13) 



rather than that of Eq. fll.2| ) the potential energy is given by Eq. ([TTT]) with n = \d s u\, and 
u is constrained to lie in a Riemannian manifold (see [0]). The resulting Euler-Lagrange 
equations for u are just second order in space and time, and the Lagrange multipliers 
can be expressed in terms of the second fundamental form of the target manifold applied 
to first derivatives of u. Both local and global existence results have been obtained for 
various wave map problems 0, |], ||, |7|, |§]. 

We note that choices for the potential energy of a geometric nature other than that of 
Eq. ( |1 . 1| ) are possible. Let n(-,t) and b(-,t) be the standard normal and the binomial 
to the curve x(-,t), as defined by the Serret-Frenet formulas, and let A = A(s,t) be 
the 3 x 3-unitary matrix whose columns are u(s,t), n(s,t) and h(s, t) (see Eq. ( [2.2QD 
below). An alternate choice for the potential energy related to the Dirichlet form energies 
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for harmonic maps [g [ig, [TTj, |T3p 

1 /" 27r ds 

Vi = - / tr (M f (M) ^(M)) - 
4 Jo 271 

= l -j\^^t) + e\s^t))^- (i.i4) 

where k is the curvature, and 9 is the torsion of the curve. Combining this with the kinetic 
energy in Eq. ( |1.2[ ) results in adding a term of the form 

d 2 ( ^d s u x d 2 s u) + d 3 s ( ?4d s u x d 2 s u) (1.15) 



w s / 1 S \ O 



to the right hand side of Eq. (|l3|), which is 6 th order in the spatial derivatives of u. 

For a physical closed loop of wire with a small circular cross section of radius p > 0, 
a more realistic expression for the elastic potential energy is given by 



7T 



A /i /" 27r 2 ^ S f 1 



V 2 = -p 4 ( T + 9/ «(M)V + 5 / #M)tt- +«(p 4 ), d-16) 



2 14 2' L K 1 ' 2tt 2 



2tt 



(is 



2tt 



where A and are the Lame constants for a homogeneous isotropic hyperelastic material 



(see [|T4|]). (For a wire made of material of a fixed density, the kinetic energy is of order 
p 2 . The kinetic energy due to twisting is of order uJ 2 p A , where u is the angular velocity.) 
The torsion term accounts for the expense of twisting the material frame. The equation 
is derived under the "quasi-static" assumption that the material arranges itself instanta- 
neously about the central curve such that the contribution of the local twist to the elastic 
energy is as small as possible. The minimizing configuration of the material for a given 
curve is achieved by a local twist which is constant along the wire. For a dynamical elas- 
tica, that constant will in general change over time. The corresponding potential energy 
for a very thin, narrow ribbon would contain a torsion term proportional to the last term in 
Eq. fll.l4| ). Maddocks and Dichmann [JT3p, Coleman et al. [ |T6[ ] and others consider direc- 



tor theories, originated by Kirchhoff and Clebsch, in which there are further stress-strain 
relations between the tangent and two independent normal vectors (see [|T7]j). The above 
choices for the potential and kinetic energy will not be pursued here. 

Our equation is related to the Localized Induction Equation (LIE) first discussed by 
Da Rios | JT~8| ] and rediscovered by Arms and Hama [JT9|] . This is a second order equation 



for the approximate time evolution of a thin vortex filament x(s, t) (e.g. a smoke-ring) 
moving in a fluid, and is given by 

d t x(s,t) =d s xxd 2 Ms,t), (1.17) 

again, s being arclength. (The equation ignores long-distance effects of the vortex acting 
upon itself.) Differentiating this equation with respect to s, one obtains an equation for 
the tangent vector u to the vortex, 

d t u(s,t) = ux d 2 u(s,t), (1.18) 



5 



known as the Landau-Lifshitz equation for the continuum Heisenberg ferromagnet [|ZU|]. 
We will elaborate on the connections between our equations and particularly the Landau- 
Lifshitz equation in the next subsection. A convenient reference for the LIE equation can 
be found in Newton [ |H| ] . 

The equations of motion for the elastica exhibit a rich variety of solutions. Langer 
and Singer found a countably infinite number of equilibrium configurations which are 
contained in tori of revolution and represent all but a finite number of torus knots, and 
showed that, up to the symmetries of the equations, every non-planar equilibrium config- 
uration appears in their list [p2|]. They also relate these solutions to those of LIE [|23|]. 

At least in modern times, consideration of small amplitude vibrations of a rotating 
(in general extensible) ring seem to go back to Carrier's 1945 paper [p4]]. He moreover 
considered the case in which the ring was constrained, or supported at points around the 
ring. Simmonds [ ^5] ] also considered small planar vibrational modes for a nearly circular, 
extensible ring, in particular flexing modes in which extension is essentially negligible. 
His analysis provides a systematic treatment of the small amplitude approximations, and 
he shows for example that the vibrational frequencies decrease with amplitude. In a 
different direction, Coleman and Dill [ ]Z5| ] examined the infinite length planar elastica, and 
showed that the solitary waves are of the form of a single loop and that traveling waves 
are a succession of periodically spaced loops, all of which satisfy differential equations 
similar to Euler's equations for the static case. Note that their equations can include 
rotational inertia, cf. their equations (44a,b). They also find wave solutions that are 
periodic in time. Following this work, Coleman and Xu [ |27| ] numerically investigated 
solitary waves for an elastica of large length, and showed that the scattering was more 
than a simple phase shift, thereby providing compelling evidence that the elastica is not 
completely integrable. In [|16|], Coleman et al. considered the elastica moving in R 3 , 
showing existence of traveling and solitary waves exhibiting torsion so that the resulting 
curves are corkscrew-like. They include some discussion of the (finite) closed curve 
case. In their work on global existence of solutions for the planar elastica, Caflisch and 
Maddocks [[[]] also showed Liapunov stability for solutions near isolated relative minima 
of the potential V. 



1.3 Description of the proof 

A first step towards solving Eq. Q1.5D is to obtain an equation for the Lagrange multiplier A 
in terms of the solution u. In the related case of a wave map, an explicit expression 
for the Lagrange multiplier in terms of first derivatives of the solution is obtained by 
projecting the equation onto the normal of the target sphere [EH. In our problem, taking 
the inner product of Eq. (JTT5J) with u(s, i), writing |d s u(s,t)| = n(s,t), and using that 
u-<9 s u = u-<9 t u = due to the constraint yields an elliptic boundary value problem for 
the tension A(-, t) at time t, 

(-d 2 s +K 2 )X = \d t u\ 2 - u • d 4 s u . (1.19) 
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A slightly less obnoxious form of this equation (involving lower order derivatives) results 
if one uses the identity 



3|d 2 u| 2 + 4<9 s u ■ d 3 s u + u • <9 4 u = , (1.20) 

which follows by differentiating the constraint |u| 2 = 1 four times with respect to s. With 
this identity the tension equation can be written as 

(-d 2 s + k 2 ) (A + 2k 2 ) = \d t u\ 2 + 2k* - \d 2 u\ 2 . (1.21) 

We do not know how to make sense of Eq. ( |1.21| ) without requiring at least u(-, t) £ H 2 , 
<9 t u(-,t)eL 2 . 

Rewriting Eq. ( |L5| ) as a system, we want to solve 

dtn = V (122) 

d t v = -«9 4 u + «9 2 (Au), u - zz; 



where A is determined by Eq. fll.21| ). The linear part of Eq. ( |1.22D generates a strongly 



continuous semigroup on H r+2 x H r for any r £ R. Standard semilinear theory accom- 
modates a nonlinearity that defines a locally Lipschitz continuous map from this space to 
itself as a perturbation. Since u £ H r+2 , it would suffice to show that A £ H r+2 to apply 
this technique. However, we only have k £ H r+1 , hence k 2 £ H r+1 (provided r > —1/2 
ensuring that H r+1 is an algebra) while A + 2k 2 is more regular by Lemma ^3] . Thus, we 
can only expect A £ H r+1 and so Eq. fll.5| ) cannot be solved directly just by converting 
it to a Duhamel integral representation. A similar picture emerges if A is inserted into 



Eq. 



Instead, we proceed as follows. Since Eq. fll.2lD is the projection of Eq. ( |L5| ) in the 



direction of u, we combine it with the complementary projection 

— u x (u x d 2 u) = u x (u x <9 4 u) + 2(d s X)d s u — Au x (u x d 2 u) . (1.23) 

We have used that u x u = and u • d s u = 0. The last two terms on the right hand side 
are locally Lipschitz from H r+2 x H r to H r . Now the problem is that the projection of 
the fourth derivative onto the orthogonal complement of u is a complicated nonlinear op- 
eration. Eqn. Q1.23D is related to the Landau-Lifshitz equation, Eq. ( [1.18D . This becomes 



apparent if we use Eq. dTTSh to formally write a differential equation for the second time 



derivative of u. The resulting equation differs from our Eq. Q1.23| ) in the last two terms: 
for us the tension is given by the solution to an elliptic problem, hence these terms are not 
local, while for the Landau-Lifshitz case the tension is given by a local expression of u 
and its derivatives. The similarity between the equations suggests how to proceed. 

The Landau-Lifshitz equation is equivalent to a nonlinear Schrodinger equation via 
the so-called Hasimoto transformation ([^8fl, cf. [[291]), which is is a nonlinear, solution- 
dependent change of variables where all partial derivatives of a curve are expressed with 
respect to a local coordinate frame transported along the curve. In Section ^, we will use 
the Hasimoto transformation to transform the projection of the second derivative operator 
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on the left hand side and the fourth derivative operator on the right hand side of Eq. (|1 .23[) 



into linear differential operators plus a perturbation, while leaving the other terms essen- 
tially invariant (see Eq. (j277l)). Unless the motion described by Eq. ( ]L5| ) is planar, the 
Hasimoto transformation introduces a monodromy into the problem. In other words, the 
frame transported along the elastica is not periodic but rather quasiperiodic, with a gen- 
erally time-dependent rotation of phase 2nP(t) about the tangent equal to the integral of 
the torsion over the curve (see Eq. ( [2.25D ). We correct for the monodromy by performing 
an additional gauge transformation. Interestingly, the Hasimoto coordinate frame with 



the monodromy correction coincides in our problem with the natural material frame [[23]] . 



However, monodromy does not arise in the Kirchhoff-Clebsch director theories [ |15| , |16U - 
Our Hasimoto transformation allows us to rewrite Eq. ( \l.5\ ) in the form 

d 
Jt 



-Y(t) = G m Y(t) + F m (Y(t)) , (1.24) 



where the linear part Gpi t ) is a differential operator which generates a strongly continuous 
evolution operator on a Sobolev space y of periodic functions, and Fpu) is a nonlinear 
perturbation. The operator GW) depends on the on the monodromy, which is in turn 
determined by an auxiliary equation 

j t (3(t) = B(Y(t)) , (1.25) 



see Eqs. (|2~TT| ) and ( p~T] ). 



In Section 0, we provide the basic estimates for the linear evolution operator Vp gen- 



erated by Gp(t) and for the nonlinearity Fp^ ■ In Lemma p3j we prove a positive lower 
bound on the spectrum of the Schrodinger operator L K = —d 2 + k 2 appearing on the 



left hand side of Eq. fll.21| ) which may be of independent interest. We conjecture that the 



lowest eigenvalue of L K is minimal when the elastica is a circle (see [BO] for a related 
result on —d 2 — k 2 ). 



In Section we set up a contraction mapping argument for Eqs. ( |1.24| )-( |Q5| ). A 



technical complication is that the dependence of the linear evolution operator Vp gener- 
ated by Gpt t ) is only strongly continuous, not norm-continuous in y r x R. We overcome 
this difficulty by setting up the contraction mapping argument using a weaker norm. 

Our local existence proof does not incorporate the more modern space-time methods 
such as Strichartz inequalities. Versions of these of these inequalities adapted to problems 



with periodic boundary conditions [31] have been used to obtain global existence of so- 



lutions for other related wave equations, for example the Boussinesq equation on a circle, 



which is also fourth order [f32p. The methods do not immediately apply here however. 
The (spatial) derivative in the nonlinear term of the Boussinesq equation is of lower order 
than in our wave equation Eq. (Oh. Thus the Duhamel integral form of the equations of 



motion analogous to our Eqs. ( |1.22P make sense for example in H r+1 x H r for suitable r, 
unlike the situation here which necessitates the Hasimoto change of variables. But again 
this change of variables leads to the monodromy issue and an associated linear problem 
with a time-dependent generator for which we do not have suitable space-time estimates 
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(see Eq. ( |1.24[ )). Finding such estimates could be an avenue to relaxing our regularity 
assumptions for the initial data. 

For the dynamical elastica moving in three dimensions, the question of global exis- 
tence remains open. The picture which does emerge from our approach is that we can 
always integrate forward for an open interval of time, i.e. we have existence (and even 
uniqueness for planar motion), until || (u(-, t), d t u(-, t))\\ H2xL 2 becomes unbounded. But 
at this moment, the tension A(s) given by Eqs. ( |1.2ip becomes infinite at some point s (or 
at least the individual terms on the right hand side of Eq. ( |1.21| ) are not integrable so that 
it is by no means clear that at this moment A exists even as a distribution). Seemingly the 
elastica would break apart. It would be of interest to know whether indeed infinite tension 
can develop in a finite time. 

We conclude with a couple of remarks about the infinite length elastica. A lower 
bound on the spectrum of L K is by no means apparent in this case, for example L K acting 
in L 2 (R) is typically not invertible, i.e. there is an infrared divergence and the estimates 
in Lemma |33] would fail. Physically, this divergence corresponds to elastica configura- 
tions with infinite tension. A local existence proof for the infinite elastica would thus 
presumably involve more subtly defined function spaces for which the tension is finite. 
The second question concerns the role of the monodromy introduced by the Hasimoto 
transformation. Even for finite-energy solutions, it is not obvious how to remove the total 
twist by a gauge transformation. 



2 A Change of Variables 

2.1 Hasimoto transformation 

Let u(s, t) be a smooth solution of Eq. (jOb. We will express the partial derivatives of 
u in terms of a positively oriented orthonormal frame which consists of u and two other 
unit vectors which we combine to a single complex vector v. The vector v(s, t) is chosen 
so that for any fixed time t, 

d s v = -(v-<9 s ii)u, (2.1) 

that is, the real and imaginary parts of v(-, t) are moved along the curve u(-, t) by parallel 
transport on § 2 . Then 

f d s u = Re [qv] ( d t u = Re [pv] 

[ d s v = —qu \ d t v = —pu + iav 

where p(s,t) and q(s,t) are complex-valued, and a(s,t) is real-valued. Since d s d t u = 
d t d s u and d s d t v = d t d s v, we must have 

(d t - ia)q = d s p, d s a = Im [qp\. (2.3) 
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In these coordinates, we compute for the projection of Eq. ( |1.5[ ) in the direction of u, 

~\P\ 2 = d] (2\q\ 2 + A) - \q\ 2 (\q\ 2 + A) - \d s q\ 2 , (2.4) 
see Eq. ( |1.21| ). For the complementary projection, we find 

(d t -ia)p = d 2 s (-d s q) + 2d s ((\q\ 2 + X)q) - (\q\ 2 + X) d s q + Re [qd s q}q , (2.5) 



see Eq. ( J1.23D . If u is a smooth solution of Eq. (]L5|), and v, a, p, q, are defined by 
Eqs. (|2.2[), and if we set 



H = X + 2\q\ 2 , (2.6) 

where k = |g| is the curvature, then we arrive at the system 

(d t — ia)p = —d^q — 4Re [q d s q]q — ilm [q d s q]q 
+2(d s /j)q + /jd s q 

(d t - ia)q = d s p 
d s a 

I (~d 2 + |g| 2 )/i 

Here, the first equation follows from Eq. (E3h, the second and third are the consistency 



(2.7) 



Im 

\p\ 2 + |g| 4 - |<9 s g| 2 



relations in Eq. (|2_3j), and the last follows from Eqs. ( |2.4| ) and (|2.6t). 



Remark 2.1 There are many other ways to complete u to an orthonormal frame, which 
are all related by gauge transformations 



J-y(s,t)~ 



e n(s,*)p ? q 



e i^' u >v , p 

a i— > d t j(s, t) + a , /i i— > ji 



e i-y(s,t)q 



(2.8) 



with some real- valued function j(s, t). The choice in Eq. Q2.1D has the special prop- 
erty that no second derivatives appear on the right hand sides of Eq. (|2.7|). As discussed 
in connection with Eq. fll.23| ), this is a key step towards solving Eqs. (|l3[)-(|T77[), because 
expressions containing first derivatives of q can, but expressions containing second deriva- 
tives of q cannot be treated as perturbations of the third derivative operator in the first line 
ofEq. Q. 



2.2 The monodromy correction 

Even when x(-,t) is 2n -periodic, the frame (u, v) defined by Eq. ( [2.1D is in general 
quasiperiodic, 



v(s + 27r,t) = e 2wim v(s,t) , 



(2.9) 



where (5 is a real-valued function of time. The value of (3(t) is determined modulo an 
integer by the monodromy of the parallel transported frame, which is related to the curve 
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torsion, see Eq. ( |2.25| ). It follows that p, q, and a satisfy time-dependent quasiperiodic 
boundary conditions. 

Periodic boundary conditions are recovered by the gauge transformation 

v(s,t) = e ^W-^v(s,t) 

P(s,t) = e^- s ^p(s,t) 

Q(s,t) = jlP®->f>®)q{ 8i t) ' 

a(s,t) = a(s,t) + d t ((3(t) - s(3(t)) , 

where (3(t) and f3(t) are chosen so that v is periodic, and a averages to zero over a period 
of 27r. We arrive at the system of equations 

(d t -ia)P = -{d s + if3fQ-ARe[Qd s Q]Q-i\m[Qd s Q\Q 

-if3\Q\ 2 Q + 2(d af j)Q + n (d s +if3)Q 
(d t -ia)Q = (d s + ip)P (2.11) 

d s a + O t p = lm[QP] 

;-^ + |Q| 2 )/i = \P\ 3 + \Q\*-\{d a +il3)Q\ 2 - 

Here, P(s,t),Q(s,t) G C, and ct(s,t), (3(t), fi(s,t) G M satisfy 2ir -periodic boundary 
conditions, with the constraint 

r 27r d s 

J a(M) 2^ = °- (2 ' 12) 

Initial conditions are given by 

P(s,0) = P (s), Q(s,0) = Q (s), 0(0) =A>, (2-13) 

where P and Qo are periodic complex-valued functions, and @ is the monodromy at 
time t = 0. We emphasize here that it is really this system, Eqs. ( |2.1 1| )-( [2~T3| ) that we 
analyze in this paper. We will see in Section |3] that Eq. Q2.1 ip can be written in the form 
of Eqs. (|P3))-(|L231) with Y = (P, Q). 

Remark 2.2 Eqs. Q2. 1 l| )-< pTT2| ) and the periodic boundary conditions are invariant under 
the gauge transformations 

Q ' \ V ' (2.14) 

pi— i>p — fc, // i— > /i , 

when fc is an integer and s G R . 

We turn to the relation of the initial- value problem in Eqs. ( |2.1 1| )-( |2T3| ) with the 
initial-value problem in Eqs. fll.5|)-(|L7|) posed in the introduction. We say the initial 
conditions (P , Qo, Po) for Eq. ( [2.1 ID are compatible with Eq. dOl), if the linear system 

d s u = Re [Qov] 2 15 

d s v = -Q Q u - i/3 v , 

has a 27r-periodic solution forming an orthonormal frame. 
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Lemma 2.3 For r > 0, the following statements are equivalent: 

1. For each pair of initial values (uq, ui) G H r+2 x H r satisfying the condition in 



Eq. f |/.(S[ ), the initial value problem in Eqs. (\l.ty - ftT7ty has a solution (u, c^u) G 
H r+2 x H r , defined on some short time interval, which is strongly continuous in 
t and assumes the initial values at t = 0. The solution is unique and depends 
continuously on the initial values. If the initial values are smooth then the solution 
is smooth in both variables. 

2. For each triple of initial values (Pq, Qq, Pq) G H r x H r+l x M. which is compat- 
ible with Eq. < \l.ty , the initial value problem in Eqs. ( 12. 11 | j-( |Z73| j has a solution 



(P,Q,P) G H r x H r+ x R, defined on a short time interval, which is strongly 
continuous in t and assumes the initial values att = 0. The solution is unique and 
depends continuously on the initial values. If the initial values are smooth then the 
solution is smooth in both variables. 



Remark 2.4 The regularity assumption r > is needed only to ensure that the right hand 
side of the equation for p in the fourth line of Eq. ( [2.1 1[ ) makes sense as an L 1 -function. In 
the proof of the lemma, we will construct a transformation (u, d t u) i— > (P, Q, (3) between 
H r+2 x H r and H r x H r+1 x M. which is continuous and has a continuous inverse for all 
r > -1. 



Proof of Lemma |2.3| : Given initial values (u , ui) for Eq. dO| ) satisfying the con- 



dition in Eq. we determine initial values (Pq, Qq, fto) for Eq. ( [2.1 1| ) by choosing a 
complex vector Vo(0) whose real and imaginary parts complement u (0) to a positively 
oriented orthonormal basis of IR 3 , and then solving Eq. ( |2.1| ) with u(s,0) = u (s) to 
obtain v (s). The initial monodromy /3 is determined up to an additive integer by 

v (2tt) = e 2 ^°v (0), (2.16) 

and Pq and Q are given by 

P Q (s) = e-^°v (s) • Ul (s) , Q Q (s) = e -^°v ( S ) • d s u (s) . (2.17) 
This defines a continuous map from H r+2 x H r to H r x H r+1 x R for any r > 0. We have 



seen above that this transformation maps smooth solutions u of Eq. (|1 .5|)-([Lq) to smooth 
solutions (P, Q, (3) of Eq. ( [2.1 1| )-( |2TT2| ). Since smooth functions are dense in H r+2 x H r , 
the transformation can be extended continuously to all of H r+2 x H r . 

Conversely, given (P , Qo, fio), let (u , v ) be a solution of the linear differential equa- 
tion Eq. ( [2.1 5D which defines a periodic orthonormal frame, and set 

Ul = Re[P v ] . (2.18) 
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By construction, (u , ui) satisfy Eq. (|1.8|). This defines a continuous transformation from 
H r x H r+1 x R to H r+2 x H r . Given a smooth solution of Eqs. (gTT|)-(glg), we can 
define a frame (u, v) by solving 

d s u = Re [Qv] f d t u = Re [Pv] (2 ^ 

9 s v = — Qu — z/5v , ^ 9tv = — Pu — ittv 

with initial conditions u(s, 0) = Uq(s), v(s,0) = v (s). Note that the third equation 
in Eq. ( [2.1 ip ensures that the two systems in Eq. ( [2.1 9[ ) can be solved simultaneously. 
Since u and v are periodic, u(-, t) and v(-, t) are periodic for t > because the pair 
of equations on the right hand side of Eq. ( |2.19[ ) preserves periodicity by the periodicity 
of P and a. The function u obtained in this way is smooth and solves Eqs. (|1.5|)-(|L7|). 
As above, the transformation can be extended continuously from the subset of smooth 
functions to all of H r x H r+1 x R . 



2.3 The standard normal frame 

It is instructive to express v, a, P, and Q in terms of standard normal coordinates, which 
are defined for smooth curves at any point where the curvature does not vanish. Assume 
u(s, t) describes the unit tangent vector of such a curve x(s,i), let n(s, t) be the unit 
normal to the curve in the direction of d s u(s, t), and set b = u x n. The vector n is called 
the principal normal, and b the binomial of the curve at s, at time t. The standard normal 
frame (u, n, b) is characterized by the Serret-Frenet differential equations 

d s u = K n 

d s (n + ib) = —Ku — i9(n + ib) 

where k(s, t) is the curvature, and 9(s, t) is the torsion of the curve. By definition, k is 
nonnegative, and 6 is real-valued. The curvature can be expressed in the various frames 
as 

K = \d s u\ = n-d s u = \q\ = \Q\ , (2.21) 

and the torsion as 

6 = U ' i9 £*f aU) =b-d s n =lm [q- 1 d s q] = Im [Q~ l d s Q] + (3 . (2.22) 
If we set 

7 (s,t) = / 6(s',t)ds' (2.23) 



(2.24) 



and define the functions v(s, t), a(s, t), p(s, t), and q{s, t) by 

v = e n (n + ib) p = e n d t u ■ (n + ih) 

a = 9(7 + n • d t b q = e n n , 
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then the differential equations in Eq. and the consistency conditions in Eq. ( |2.3| ) are 
satisfied. The monodromy (3 can be defined by 

f 2lT ds 

(3{t) = J e(s,t)—; (2.25) 

in particular, the monodromy is an intrinsic quantity associated with the curve. 

For planar curves, the Hasimoto transformation and the analysis of the transformed 
initial value problem simplify considerably. For a curve lying in the xi-x 2 -coordinate 
plane, we choose Re [v] to be the unit vector obtained by rotating u counterclockwise 
through an angle of ir/2, and Im [v] to be the unit vector in the x 3 -direction, and set v = v. 
The partial derivatives of u and v satisfy Eq. ( |2T2"| ) with a = and p, q real-valued. The 
equations on the left hand side of Eq. ( |2.2| ) agree with the planar Serret-Frenet equations, 
and q(s,t) is the signed curvature of the curve x(s,t). In this case, no monodromy 
correction is required and (P, Q) = (p, q) satisfy Eqs. Q2.1 ID with a — [3 — 0. The 
resulting system can be written as a semilinear equation 

jY(t) = GY(t) + F(Y(t)) (2.26) 

in a suitable Sobolev space of periodic functions, which can be solved by a standard fixed 
point argument (see Corollary pi~2| ). 



3 Estimates 

In the previous section, we have changed variables in the equations for or the dynamical 
Euler's elastica, and transformed Eqs. (|T31)-(|T77l) into the equivalent initial value prob- 
lem given by Eqs. ( |2.1 1| )-( |2T3| ). To see that the resulting equations have the form of 
Eq. ( |1.24{ )-( [TT23| ), set Y = (P, Q), and define the linear operator by 

Gm \Q) = \{d. + W)) o ){q)> (3 - 1} 

where (3 is is a real-valued function of time. Suppressing the time-dependence in the 
notation, we write the nonlinearity in Eq. ( |1.24| ) as a sum of three terms: 



F,(Y) = i«g) + ( 2Mg+ * +iW ) 



+ 



-4Re [Qd s Q]Q- Urn [Qd s Q]Q - i(3\Q\ 2 Q 




FW{Y) + F® (Y) + F®(Y), (3.2) 



where a is determined by 

"2tt . 



d s a = Im [PQ] - Im [P(s)Q(s)] ^ , J a(s) ^ = , (3.3) 
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and \i solves the elliptic boundary value problem 

(|Q| 2 -<9 2 )/i = |P| 2 +|Q| 4 -|(«9 S + ^)Q| 2 (3.4) 
with periodic boundary conditions. The nonlinearity in Eq. ( |1.25[ ) is given by 

BOO = I Im [Q(s)P(s)} - . (3.5) 

3.1 The linear part 

We begin the analysis of the system in Eqs. ( [2.1 1| )-( [2T2"| ) by solving the linear equation 



Jt (q) ~ G ^ (q 



G m ( U ) (3-6) 



with a given Lipschitz continuous function (3 defined on [0, T]. The fundamental solution 
of Eq. will be denoted by Vp(t, t ). 

It is natural to consider Eq. ( P.6Q in the Fourier series representation 



±(P(n,t)\ ( *(n + /3(t)) 3 \ (P(n,t)\ 

dt\Q(n,t)J " \i(n + /3(t)) J \Q(n,t) J ' 

where it decouples into a sequence of linear ordinary differential equations on C 2 . Let 
Gp(n, t) be the 2x2 matrix appearing on the right hand side of Eq. (|3.7|). The fundamental 
solution Vp{n, t, t ) of Eq. ( |3~77| ) is given by the time-ordered exponential of Gp{n, t). 

Denote by H r the Sobolev space of 27r-periodic complex-valued functions (or distri- 
butions, when r < 0) having r fractional derivatives, with norm 



(l + n 2 Y\f(n)\ 2 . (3.8) 



Let y r be the space of 27r-periodic functions Y = (P, Q) in W x H r+1 , with norm 

WYfyr = WPWht + WQWW- (3-9) 
For a vector (a, b) in C 2 , we define its n-norm by 

l) = |a| 2 + ^H 2 |6| 2 , (3.10) 

/ n 

where w(n) = y/1 + n 2 . With this notation, Eq. becomes 

oo 

WYfyr = J2 w ( n ) 2r \\Y(n)\\ 2 n , (3.11) 

Tl— — OO 

where Y{n) = (P(n), Q(n)) G C 2 is the Fourier transform of Y. 
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Lemma 3.1 Suppose that 3 is a real-valued Lip <schitz- continuous function on [0, T] with 
Lipschitz constant r], and that | (3(0) | < 1. There exists an increasing continuous function 
C of two variables with sup^ C(i], 0) < oo such that for < to < t < T and any value 
ofr E E, 

\\Vp{t,t )\\yr<C(ri,T). (3.12) 
Moreover, Vp{t, t ) is strongly continuous in y r with respect to t and t . 

Proof: We first show that 

\\V (n,t,to)\\n<C(Ti,T) (3.13) 

uniformly in n, with C as in the statement of the lemma. The idea is to rewrite Eq. ( p.7| ) 
as 



d 



P(n,t) 



dt \(n + (3)Q(n,t) 



(3.14) 



1 



P(n,t) 



+ 



dp 



i(n + (3) ^ ; V( n + ^)Q(n,t); ' dt \0(v.t) 
Let Up(n, t, t ) be the unitary 2x2 matrix defined by 

U n ((3,t,t ) = expL j (n + P(t')) 2 dt' J 



(3.15) 



By the Duhamel integral formula, Eq. ( [3 . 1 4[ ) is equivalent to 



P(n,t) 
{n + (3)Q(n,t) 



* A , i\ dB/dt 
U p {n,t,f 



to 







+ Up(n,t,t 



(n + 3) \(n + 3)Q(n,t') 

P(n,t ) 
n + 8)Q(n,t 



dt' 



(3.16) 



provided that \n\ is sufficiently large so that n + (3(t) does not vanish anywhere on [0, T]. 
Since Up(n, t, t ) is unitary, and |/3(t)| < 1 + i]T for < t < T by assumption, we can 
estimate for |n| > 2(1 + r/T), 



P(n,t) 
[n + p)Q(n,t) 



< 



2r/ 



\n\ 



+ 





n + 3)Q(n,t' 

P(n,t Q ) 
{n + 3)Q{n,t 



dt' 



(3.17) 



Applying Gronwall's inequality, and using the fact that the left hand side is equivalent to 
the n-norm, we arrive at 



V (n,t,t o ) <6e 2 ^/l"l, (M>1 + 77T). 



(3.18) 
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For any value of n, we can bound the n-norm of Gp(n, t) of Eq. ( |?77| ) on C 2 by 

Gp(n,t) 

n 

< max I (|n| + 1 + r)T)w(n) 



sup 



i{n + f3{t)f\ fa 

i{n + (3{t)) ) \b 

(\n\ + l + r]T) 3 



w(n\ 



Estimating the right hand side and applying Gronwall's inequality gives 



Vp(n,t,t ) 



< e 



2(\n\+l+r 1 T) 2 (l+r l T)T 



(3.19) 



(3.20) 



Combining Eqs. Q3.18D and ( |3.20D implies Eq. (claim :Vbetabound-n) with the value of the 
constant given by C(r], T) = Qe 8 ^ 1+riT ^ T . This proves the claim in Eq. ( |3. 12| ). 

Clearly, each Fourier coefficient Vp(n, t, t )Y(n) depends continuously on t and t . 
Since Eq. ( |3.13[ ) implies a uniform tail estimate on ||V/j(n, t, to)Y(n)\\, it follows that 
Vp(t, to) is strongly continuous in both time variables. 



We also need to bound the dependence of Vp on (5. 

Lemma 3.2 Assume that for some T > 0, the functions (3\ and fii are Lipschitz contin- 
uous on [0, T] with Lipschitz constant r\, and that |/3i(0)|, |/?2(0)| < 1. There exists an 
increasing continuous function C of two variables such that 



^ a (t,to)-^(t,to)||vr^-i < C(rj,T) / \^ 2 (t') - ^{t')\dt' 



(3.21) 



Moreover, Vp is strongly continuous in [3 with respect to the y -topology in the sense that 
for every 7 6 7 and any given e > there exists 6 = 5(e, i], T, Y) such that 



implies 



sup |/? 2 (t) < 5 

0<t<T 



sup \\(Vp 2 (t,t Q )-Vp 1 (t,t ))Y\\ < e. 

0<t <t<T 



(3.22) 



(3.23) 



Proof: We first show that for < t < t < T, 

Vp^t^-Vp^nXto) < C( V ,T)w(n) I \p 2 {t') - f3 1 {t')\dt' 



(3.24) 



to 



uniformly in n, which clearly implies Eq. ( p.21| ). To see Eq. ( |3.24[ ), we write 



Vp 2 {n,t,t ) -Vp^n^ta) = / Vp^n^AGin^Vp^t' ,t )dt' 



(3.25) 



to 
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where 



i(n+f3 2 (t)f -i(n+A(t)) 3 



AG(n,t) = 

Using that \/3i(t)\ < 1 + 77T by assumption, we estimate 



AG(n,t) < \fo(t) - lh(t)\ 



3(|n| + 1 + r]Tf 

1 



< max ^ 3 (|n|+1 ^ r)2 , W (n) - /3i(t)\ 



win) 



< 3(1 +r ] T) 2 w(n)\P2(t)-p 1 (t)\. 



(3.26) 



(3.27) 



Inserting Eq. ( |3.27[ ) and the bound in Eq. ( [3.13D of Lemma [TI] into Eq. ( |3.25t ) yields 
Eq. (gggj), with C(r/,T) = c(r?,T) 2 3(l + r]T) 2 , where c(r/,T) is the constant from 
Lemma |3T|. This proves the claim in Eq. ( |3.21| ). 

The strong continuity follows as in the proof of Lemma |3~T| by combining Eq. ( [3.24D 
with a uniform tail estimate obtained from Eq. ( P-12[ ). 



3.2 The resolvent for — dl + k 



We provide a lower bound for the spectrum of the Schrodinger operator L K = — d 2 + k 2 
acting in L 2 [0, 2tt], with periodic boundary conditions. 



Lemma 3.3 The operator L K has least eigenvalue eo(/«) satisfying, 

eo(«) > 1/4 . 

/n particular, L K is invertible in L 2 , and IIL" 1 ^ < 4. 



(3.28) 



Remark 3.4 It is natural to conjecture that inf 6q(k) is actually attained for x(s) a circle 
where k(s) = 1 and e (n) = 1. An analogous result due to Harrell and Loss says that the 



second-lowest eigenvalue of — d s — k is maximal for a circle [pOQ. 



Proof of Lemma [O]: We have that 



inf 



inf + 

{i/):||V'||2=l} 



> inf '||9 S *||2 



(3.29) 
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where the vector function ^ = ^u, and in the last line, the infimum is to be taken over 
all normalized ^ such that 

r \&( s ) 

u(s) ds = / —Hi ds = 0. (3.30) 
J 1*001 

If these constraint integrals are zero, then each component ipi of ^ must vanish at some 
point i = 1,2, 3, i.e., each component must satisfy a Dirichlet condition and so 
II dsV'i I ! > IllV'illl' f° r eacn * so tnat tne infimum in the last line of Eq. ( |3.29j ) is bounded 
below by 



mf y^iiaWill! > inf -Y UiWi = -• ( 3 - 31 ) 



This completes the proof of the eigenvalue estimate. 

■ 

Next, we provide bounds on the resolvent for L K , considered as mapping H r to H r+2 , 
assuming the spectral bound in Lemma |3~3| and additional bounds on the norm of k. The 
inequalities given are by no means optimal, but they are adequate for our purposes; the 
estimates are in the spirit of Bessel kernel estimates (see [03|]). 



Lemma 3.5 If k G H r+1 for some r > 0, then L~ x defines a bounded linear operator 
from H r ~ 1 to H r+1 . More precisely, there exists a constant C\ = C\{r) such that if p 
solves 

(-d 2 + K 2 )p = f (3.32) 

with f G H r ~ l for some K with \\K\\ Hr+1 < R, then p satisfies 

\\p\\ H r + i < C^l + RY-WfWiir-i, (3.33) 

where v = v(r) is the smallest the smallest integer at least as large as (r + 1)/ 2 for r > 1, 
and u(r) = 2 for < r < 1. The Fourier coefficients of p are bounded by 

\p(n)\ < C 2 w{n)- 2 (^{n^+winy^il + R 2 )" -\\f\\ Hr ^ . (3.34) 

Proof: We compare L K with the operator —d 2 + 1, using two forms of the resolvent 
identity: 

L- 1 = (-d 2 s +iy l -(-d 2 s + l)-\K 2 -l)L- 1 (3.35) 
= (-^ + l)- 1 -^V-l)(-^ + l)- 1 - (3-36) 

Clearly, 

||(-^ + l)-V||^, +1 = 11/11^ (3.37) 
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for any r' E R. Let f E L . The first line of the resolvent identity shows that 

IfcVIU < ii/ii^ + IK^-i^/L 

< Cl (l + i? 2 )||/|| i2 (3.38) 



for some constant cj. We have used Eq. (3.37) in the first line. In the second line, we have 



used Lemma p.8| , which is proved below, and the spectral bound of Lemma |3_3|. This 
proves the claim in the case r = 1 . More generally, the first line of the resolvent identity 
shows that for r' < r, 

||^'«; 1 || j ffr'+l_ >jff r'+3 — C 2(l + H-^K^I^r'-l^^r' + l (3.39) 

for some constant C2 = C2(r, r'). Iterating this estimate we obtain the claim for r > 1. 
The case < r < 1 follows from the second line of the resolvent identity, which implies 
that for r' < r 

\\ L ^\\ H r'-^Hr'^ < ^(l + R^WL- 1 ^,^^ . (3.40) 

■ 

We will also need estimates on how /J varies as the curvature k varies. 
Lemma 3.6 Suppose that /ii and \i2 solve the equations 

L Kl m = (-<9 S 2 + = / 

L K2 fi2 = (~d 2 s + kI)/j, 2 = f , (3.41) 

where Ki,k 2 E H r+1 . Suppose that R > max {||fti||#r+i , || K 2||//r+i} • For r > there 
exists a constant C = C(r) such that 

||A*2 - A*i||fl-r+i < CR (l + R 2 ) 2v ||« 2 - K\\\ H r ll/ll^-i ■ (3-42) 
Here, v = u(r) is the exponent from Lemma 
From the resolvent identity 

/la -Mi = -L-l{K\-K\)L-lf , (3.43) 
we obtain for r > with the help of Lemma pT8j , which is proved below, 



||M2-Ml|| r+ 1 < ^H^Lr^+l ||«2 -Kl\\ H r H^lljyr-l^, ll/lltf-l -0 M ) 

The claim now follows from Lemma 
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3.3 The nonlinearity 



Lemma 3.7 For Y = (P, Q) G y r and j3 G R, to i*]g Z?<? defined by Eq. ( p72| J, and to 
-B /?e defined by Eq. There exist increasing continuous functions C\, C 2 , and C3 0/ 

variables, which also depend on r, such that the following estimates hold: 



1. If r > 0, then for any Y Ey with \\Y\\y r < R, and any f3 G R vwY/z |/3| < 6, 

|B(Y)| < i? 2 . 



(3.45) 



2. Ifr>0, then for any pair y\,Y 2 G y r with ||li||yr < -R and any /3 l5 /3 2 G R w/f/z 

IAI<&, 



iF^^-FftCn)!!^ < C 2 (b,R)(\\Y 2 -Y 1 \\ yr + \p 2 -p 1 \) 
\B(Y 2 )-B{Y X )\ < R\\Y 2 -Yi\\yr . 



(3.46) 



3. Ifr > 1/2, then for any pair Y\,Y 2 G y r with \\Yi\\y r < Rand any j3i, j3 2 G Rwzf/z 
IAI < b, 



iF^-F^Wy^ < C 3 (6,i2)(||y 2 -yi|| 3 ^ 1 + |^-/3i|) 

|s(y 2 )- J B(y 1 )| < i2||y 2 -nily-i ■ 



(3.47) 

The following lemma will be useful in estimating the various terms in the nonlinearity. 



Lemma 3.8 (Leibnitz rule for if r -norms) There exist constants C\ - C5 which depend 
only on r, so that each of the following inequalities holds whenever the right hand side is 
finite: 



(3.48) 



Proof: If r is an integer, then the first two inequalities follow immediately from the 
Leibnitz rule and the fact that H 1 C L°° in one space dimension. In general, we use that, 
for r > 0, w(n) r < Ci(w(k) r + w(n — k) r ) with some constant c\ = ci(r) to get 



\\fg\\ H r 


< 




\f\\ H r 


\\g 




(r > 1) 


\\fg\\ H r 


< 


c 2 


\f\\ H r 


\\g 


|ffr+l 


(r > 0) 


ll/^llffr-l 


< 


c 3 


\f\\ H r 


\\g 


\H r 


(r > 0) 


II/^IIjt-i 


< 


c 4 


\f\\ H r- 


-1 1 


g\\H r + i 


(r > 1/2) 


\\fg\\ H r-* 


< 


c 5 


\f\\ H r- 


-1 1 


g\\H r 


(r > 1/2) 



w(nY\Tg(n)\ < w(nyJ2\f(k)\-\g(n-k)\ 

k 

ci[J2 w ( k y\M-\g(n-k)\ + J2\Hk)\- 
V k k 



(3.49) 



w(n — k) r \g{n — k)\ 
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If either f,gE H r with r > 1, or / 6 H r , g E H r+1 , each of the two sums on the right 
hand side is the Fourier transform of the product of an L 2 -function with an if ^function, 
and hence in £ 2 . This proves the first two inequalities. 

For r > 1, the third inequality follows from the first. For < r < 1, we use that 

w(n) < y/2w(k)w(n — k) to obtain 

w{n) r ~ l \Tg{n)\ < c 2 w{n)- 1 J2 w (k) r \Kk)\-w{n-k) r \g{n-k)\, (3.50) 

k 

where the constant c 2 depends only on r. The sum is the Fourier transform of the product 
of two L 2 -functions, and hence in £°°. Since w^ 1 £ £ 2 , it follows that w(n) r ~ l fg e £ 2 , 
which implies the third inequality. 

For r > 1, the last two inequalities follow from the second and third with r replaced 
by r — 1. For 1/2 < r < 1, we use that \^2w(n) > w(k)/w(n — k) and proceed as in the 
proof of the third inequality. 



Proof of Lemma |3T7] The estimates for F (3) , and B follow by repeated applica- 
tions of Lemma [3T8| . We focus on the terms involving /i. 

Consider the first claim, Eq. fl3.45| ) for r > 0. We want to apply Lemma |3]5| with 
K = \Q\ e H r+1 and / = \P\ 2 + |Q| 4 - \(d s + i(3)Q\ 2 . By the third inequality of 



Lemma |3i| we have for r > 0, 

WfWm-i < Cl (l + \(3\ 2 )R 2 . (3.51) 

By Lemma |3~3| , 

< c 2 (l + W)R\l + R 2 Y , (3.52) 
and so by the second inequality of Lemma OL 



\F ( *\Y)\\yr = \\2(d af j.)Q + n(d s + i/3)Q\\ Hr 

< c 3 (l + |/3|)||//||^ +1 ||g||^ +1 (3.53) 

< Ci (l + \P\ 3 )R\1 + R 2 Y , 



where v = u(r) is the exponent from Lemma |3~5|. This shows the bound in Eq. (|3.45|). 



The proof of the second claim, Eq. ( J3.46[ ), is almost the same. 



To see the third claim, Eq. ( p.47[ ), let = (Pi, Qj), «^ = \Qi\, and denote by fa the 
right hand side of Eq. ( pT4[ ) corresponding to Y { (i = 1, 2). By Lemma |3T8| , we have for 

r > 1/2, 

< c 1 (l + b 2 )R\\Y 1 -Y 2 \\ y ^+c 2 bR 2 \/3 2 -p 1 \ (3.54) 

< c 3 (l + b 2 )(l + R 2 ) (UYa-nil^i + \(32-(3i\) 
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with suitable constants C1-C3 By Lemmas R3I and KR this implies 



\V>2-V>l\\ H r < \\L K2 (f2-f 



(L-l - L-l) h 



< c 4 {l + b 2 ){l + R 2 Y\\f 2 -h\\ Hr 

+c 5 R(l + R 2 ) 2l/ \\k 2 - Kx\\ Hr 

< cb(1 + b 2 )(l + i? 2 ) 2 ^ 3 / 2 (||F 2 - Y l \\ yr _ 1 + \(3 2 - Pi 



(3.55) 



with suitable constants c^-cq. Inserting this and Eq. ( J3.54[ ) into and using the fourth 
inequality of Lemma OL we see that satisfies the bound in Eq. (|3.47|). ■ 



4 Proof of Theorem [tt] 

4.1 Existence of solutions to Eqs. ( |2mD -( |27r3D 

In the notation of Eqs. (|3TT|)-(|3T6|), the Duhamel formula for Eqs. ( |2. 1 1[ )-( |27T3| ) is given by 



Y(t) = Vp(t,0)Y o + f*V p (t,t')F m (Y(t'))dt =: Fp{Y){t) 
m = (3 + JoB(Y(t'))dt' =: B(Y)(t), 

where Y = (P , Qo) £ y r and f3 G M. are given initial values. We begin by solving the 
first equation in Eq. Q4.ip for a fixed function (3. 

Lemma 4.1 Let (3 be a Lipschitz continuous real-valued function on M + with Lipschitz 
constant r\ and |/3(0)| < 1, and to lo G y for some r > 0. TTzere existe a number 
R < 00, which depends on \\Y \\yr, and a time T > which depends on 77 and R such 
that the fixed point equation 

Y = JFg(F) (4.2) 
/ias a unique solution Y£ on [0, T] m 3^ r . 77ze solution satisfies the uniform bound 

sup ||r;(t)|L r (4.3) 

0<t<T ^ 

an J depends continuously on the initial value Y with respect to the y r -norm. It is also 
strongly continuous in [3 with respect to the y r -norm, uniformly on [0, T], in the sense 
that for any e there exists 5 = 6(e, (3, T, y Q ) such that 

sup \{3'{t)-(3{t)\ < 5 (4.4) 

0<t<T 

implies that the corresponding solutions YS and YZ, satisfy 

sup \\Y*(t)-Y;(t)\\ yr < e. (4.5) 

0<t<T 
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Proof: We consider Tp as a map from the space 

V R = \y : [0,7] h-> y r I Y continuous, sup \\Y(t)\\y r < r\ (4.6) 

I 0<t<T J 

with norm 

|||y||| r = sup ||y(t)||yr (4.7) 

0<t<T 

into the space of real- valued continuous functions on [0, 7} with values in y r . The values 
of 7 and R will be chosen below. 

By Lemma |3J| , there exists an increasing continuous function of two variables C\ 
with sup^ Ci(r), 0) < oo such that for any Z £ y r , we have 

sup ||^(t,t )^||yr <Ci (77^)11^11 (4.8) 

0<t <t<T 



By Eq. < J3.46| ) of Lemma [577] , there exists an increasing continuous function C 2 of two 



variables such that for any two continuous functions Yi, Y 2 on [0, 7] with values in y r 
which are bounded uniformly by R, 

\\F m (Y 2 (t)) - F^Y^t))^,, < C 2 (|/3(t)|,it;)||y 2 (0-yi(t)||^ (4.9) 

for all < t < 7. Combining the above two estimates and using that \/3(t)\ < 1 + r/7, 
we see that 

sup \\F (Y 2 )(t) - FpiYj^Wyr 

0<t<T 

< [ C 1 ( V ,7)\\F(Y 2 (t'))-F(Y 1 (t'))\\ yr dt' (4.10) 
Jo 

< 7C 1 {r ] ,7)C 2 (l + r ] 7,R) sup \\Y 2 (t) - Y 1 (t)\\ yr . 

0<t<T 



Note that by Lemma |3J] and Eq. ( |3.46D of Lemma |377| , the function J-p(Y) * s a g am a 



continuous function of t. 
Fix 

R > 4supC 1 (r/,0) ||y |L r , (4.11) 

r?GR 

and choose 7 small enough so that 

d^T) < 2C 1 (t ?) 0), 7C 1 (r,,7)C 2 (l+r,7,R) < i. (4.12) 
Then JF^ is a contraction with Lipschitz constant 1/2 on Dr. Since, for y £ Pr, 

sup ||^(Y)(t)b- < sup \\v f3 (t,o)Y \\y r + sup ||^(y)-^(o)b, 

0<t<T 0<t<T 0<t<T 

< R, (4.13) 
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we see that Tp maps Vr into itself. By the contraction mapping principle, Tp has a 
unique fixed point in Vr, which we denote by Y*. 

The function Tp is clearly continuous in the initial value Y with respect to the norm 
on Vr. By Lemma [^2] and Eq. ( p.46[ ) of Lemma [J77| , it is also strongly continuous in (3, 



in the sense that for every e > there exists S = 5(e, rj, T, Y) such that 



implies 



sup \/3'(t)-p(t)\ < 5 (4.14) 

0<t<T 



sup \\Fp>(y) - F p (X)\\yr < e, (4.15) 

0<i<T 



where Y is any continuous function on [0, T] with values in y r , and (3 is Lipschitz con- 
tinuous with |/?(0)| < 1 and Lipschitz constant r\. By the uniform contraction principle 



(see [p4|], Theorem 2.2), the fixed point Y* inherits the claimed continuity properties from 
Tb- The modulus of continuity depends on (3 and on Y through the dependence of Yg* on 
these parameters. 



Corollary 4.2 Let Y = (P , Q ) G y r for some r > 0, where P and Q are real-valued, 
and let j3 Q = 0. There exists a time T > 0, which depends on ||Yo||y> sucn tnat the fixed 



point equations in Eq. ffO]) have a unique solution of the form (Y*, 0) on [0, T], where 
Y* = (P*, Q*) G y r are real-valued. The solution depends continuously on the initial 
values in y r . 



Proof: The claim follows immediately from LemmaO with (3 = 0. 



Theorem 4.3 Let r > 0, and let initial values Yq G y r and (3q G R be given. There 
exists a time T > 0, which depends on HYollyr |/5o|» ^"c^ that the pair of fixed point 



equations in Eq. (4_A_) has a solution in y r x R on [0, T] which is strongly continuous in 



t and assumes the initial values at t = 0. 



Proof: We may assume by the gauge invariance given in Eq. ( |2.14[ ) that \(3 \ < 1. Let 
R = R(\\ Y \\yr) be the constant appearing in the statement of Lemma (p~Tl), and consider 
the space 

C R = {(3: [0,T] ^ R | |/3(0)| < 1, Lip(/3) < P 2 } (4.16) 

with the norm of uniform convergence. Then Cr is a compact convex subset of the space 
of continuous real- valued functions on [0, T]. 

By Lemma |4TT|, there exists a time T = T(R) > such that for (3 G C R , the fixed 
point equation in Eqs. ( pT2| ) has a unique solution Y^* which is continuous in t and satisfies 

sup \\Y* p {t)\\yr < R. (4.17) 

0<t<T 
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By Eq. (gg5J) of Lemma we have that 

|B(r;)(i)-B(y;)(t )| < f B{Y;){t')dt> < R 2 \t-t \, (4.i8) 

that is, -8(1^*) is again Lipschitz continuous with Lipschitz constant R 2 , and hence in 
Cr. Since the map [3 t— > B(Y^) is continuous on V R by Lemma [TTJand Eq. ( P-46| ) of 
Lemma pTTj Schauder's theorem implies that it has a fixed point in Cr, which we denote 
by P*. By construction, the pair (YL , /?*) solves the fixed point equations in Eq. (f3~T|). 



Remark 4.4 The dependence of R and T on the initial value /3 is due to the fact that the 
gauge transformation in Eq. ( |2.14| ) may change the norm of the initial values Y . Let k be 
the integer closest to /3 , and perform the gauge transformation in Eq. ( [2. 140 with s = 0. 
Then \(3 - k\ < 1 /2, and for r > 0, 



iks\/ I 



< (l + A; 2 )^ 1 )/ 2 !!^ 



(4.19) 



4.2 Well-posedness 



In order to exploit the equivalence of Eqs. ( |2.1 1| )-( |2~T3| ) with the original initial value 
problem in Eqs. (|1.5|)-(p77|) that was established in Lemma we need to prove that the 
solutions of Eqs. 02. 1 l[ )-( pTT3| ) are unique, and depend continuously on the initial values 
in some reasonable norm. The difficulty is that the linear operator Vp is only strongly 
continuous, not norm-continuous, with respect to the parameter /3 (see Lemma |3?2|). We 
deal with this problem by using a weaker norm to estimate T$ and B. 

Theorem 4.5 Let r > 1/2. For given initial values Yq g y r and (3q G IR, there exists a 
time T > which depends on || lolly and on \/3q\ such that the pair fixed point equations 
in Eq. ( \4.1\ ) has a unique solution in y r x R on [0, T] which is strongly continuous in 
t. The solution depends continuously on the initial values with respect to the norm on 

y r x R 



Proof: By the gauge invariance in Eq. ( |2.14[ ), we may assume that |/5 | < 1- For R to 



be chosen below, we consider the right hand side of Eq. (|4JJ) as a map on the space 



V R xC R = < (Y, (3) : [0, T] i-> y r 



Y continuous, 

SUp <t<T ll^OOb- < R i 



Lip(/3) < R 2 



with the norm 



< 1 



(4.20) 



\\\(Y,f3)\\\ = sup ||Y(f)||yr-i+ sup \/3(t)\. (4.21) 

0<t<T 0<i<T 
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Note that V R xC R is complete with respect to the 1 1 1 • 1 1 |-norm by the convexity of the y r - 
norm and the fact that the uniform limit of Lipschitz continuous functions with a given 
constant is again Lipschitz continuous with that constant. 

We first bound the 3^ r -norm of Fp(Y) and the Lipschitz constant of B(Y). By 
Lemma |3~T| , there exists an increasing continuous function C\ of two variables with 
sup^ Ci(r], 0) < oo such that for any Lipschitz continuous function (3 and any Z G y r , 



sup \\V p (t,t )Z\\ Y „ < C x {Uv{f3),T)\\Z\\yr 



(4.22) 



0<to<t<T 



By Eq. ( |3.45D of Lemma [377] , there exists an increasing continuous function C 2 of two 
variables such that for any (Y, (3) G V R xC R , 



sup_||F^ ) (y(0)||^ < C 2 (\/3(t)\,R). 



0<t<T 



Combining Eqs. (ggg)-(ggg), we see that for (Y, (3) G V R x C R , 

sup \\Fp{Y,(3){t)\\ yr < C^T) ||y ||yr 

0<t<T 



(4.23) 



(4.24) 



+TC 1 (R 2 , T)C 2 (1 + R 2 T, R) 



By Eq. ( |3.45t ) of Lemma |3?7|, we also have 



Lip(£(Y)) 



sup 

0<t <t<T 



\B(Y)(t) - B(Y)(t )\ 



t-tn 



< R 2 . 



(4.25) 



Next we bound the Lipschitz constant of (J 7 ^, B) on X. Let (Yi, (3i) and (Y 2 , (3 2 ) be in 
T^r x Cij. By Lemma [?72|, there exists an increasing continuous function C 3 of two vari- 
ables such that for any two real- valued functions (3i, (3 2 which are Lipschitz continuous 
with Lipschitz constant R 2 , 

\\Vp,{t,t Q )Z-V Pl {t,t Q )Z\\ Y ^ < TC 3 (R 2 ,T) sup \[3 2 (t')-f3 1 (t')\\\Z\\ yr (4.26) 

0<t'<T 



for all < t < t < T. By Eq. ( |3.47[ ) of Lemma |3.7|, there exists an increasing continuous 
function C 4 of two variables such that 



W^))IU-! < C i (\f3\,R)(\\Y 2 (t)-Y 1 

+\fo(t)-Pi(t)\ 



(4.27) 



for all < t < T. Combining Eqs. (^2"2>(F2"31) with Eqs. (fiJty-tfTTfy, and using that 

\P(t) | < 1 + R 2 T, we see that for (Y, (3) G V R x C R , 



sup ||^> 2 (Y 2 , (3 2 )(t) - ^(Y, A 

0<t<T 



(4.28) 



< TC 5 (R,T, \\Y \\ yr )\\\(Y 2 ,(3 2 ) - {Y 1 ,(3 1 ) 
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with some increasing continuous function C5 of two variables. By Eq. ( ]3.47| ) of Lemma |3.7|, 

sup |£(Y 2 )(t)-£(Yx)(t)| < TRWYz-Y^y^ . (4.29) 

0<t<T 

In summary, Eqs. (g^)-(g^ show that 

\\\(Fp 2 (Y 2 ),B(Y 2 )) - (Fx), B{Y 1 ))\\\ (4.30) 

< TC 6 (R,T, \\Y \\yr)\\\(Y 2 ,p 2 ) - (Y^IW 

with a suitable increasing function C§. 
Choose 

R > 4supC 1 (r ? ,0) ||y ||yr , (4.31) 
v 

and T small enough such that 

d(R 2 ,T) < 2Cx(R 2 ,0), 
TC 1 (R 2 ,T)C 2 {1 + R 2 T,R) < R , (4.32) 
TC 6 (R,T) < \. 

Then Eqs. (g^)-(g^ show that (F p , B) maps V R x into itself, and Eq. (JO0|) shows 
that (Tp,B) is a contraction with Lipschitz constant 1/2. By the contraction mapping 
theorem, Eq. ( |4~T| ) has a unique solution in T>r x Cr . By the uniform contraction principle 
and the continuity properties of Va(£, t') and Fp, this solution is continuous with respect 
to the initial data (y , (3 ) in the sense that for every e > there exists 5 = 5(e, T, (3 , Y ) 
such that 

\f3 r -f3 \ + \\Y{-Y \\yr-i < 5 (4.33) 
implies that the corresponding solutions (Y, (3) and (Y', 13') of Eq. ( }4.1D satisfy 

|||y' -y||| < £• (4.34) 

It remains to prove the strong continuity of the solution with respect to Y and /5 in 
the natural norm on y r x R. Note that Eq. ( |4.34[ ) implies in particular that 

sup \(3'(t) -P(t)\ < e . (4.35) 

0<t<T 

Fix initial values (y , (3q) in y r x R with |/3 | < 1, let e > be given, and suppose that 

|^-/3o| + ||yo-y ||^ < ^ (4.36) 

Let (y, /3) G P/? x Cr be the solution of Eq. ( pTTj ) with these initial values, defined on 
some interval strictly containing [0, T], and let (Y', (5') be the solution to Eq. Q4.1D with 
initial values (yj, where Y G 3^, and |/3 | < 1. If Yq' is sufficiently close to Y" , then 
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we may assume that (Y'ft') is defined on [0, T]. By Lemma |Q|, Y is the unique solution 
of the first fixed point equation in Eq. ( [4.1[ ) with (3 fixed and the given initial value Y . Let 
Z be the unique solution of the fixed point equation 

Z(t) = V p (t,0)Y£+ [ V p (t,t')F m (Z(t'))dt' (4.37) 
Jo 

with the same function (3. Clearly, 

SUp \\Y'(t)-Y(t)\\yr < SUp (\\Y'(t)-Z{t)\\yr + \\Z{t)-Y(t)\\yr) . (4.38) 

0<t<T 0<t<T v ' 



By the continuity statement in Eq. ( |4.35| ), we may assume that sup 0<t<r \(3'(t) — (3{t)\ is 



as small as we please. Since Y' and Z solve the fixed point equation in the first line of 
Eq. ( [4. lj) with the same initial value Yq but different functions j3' and f3, and since Z and 
Y solve the equation with the same function (3 but different initial values Yq and Y , the 
continuity statements of Lemma pTT] imply the claim. The modulus of continuity depends 
on p and Y through the dependence of (Y, (3) on these parameters. 



Proof of Theorem [T7TJ By Lemma the well-posedness of Eqs. (|T3t)-([T77l) asserted 



in Theorem 4.5 for the three-dimensional case and in Corollary 073 for the planar case 



implies the claims of Theorem [TTT] . 
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